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In  the  paper , a self-consistent  iterative  kinetic  model  of  stratified  low-pressure  argon  glow 
discharge  in  plane  geometry  is  presented.  It  is  shown  that  striations  moving  from  the  anode  to  the 
cathode  side  of  the  positive  column  can  be  obtained  from  the  consideration  of  non-local 
Boltzmann  equation  for  electron  distribution  function  and  non-stationary  continuity  equation  for 
ions  coupled  with  Poisson  equation  for  electric  field. 


1.  Introduction 

It  is  known  that  the  positive  column  of  a low-pressure 
glow  discharge  in  inert  gases  often  exists  in  a stratified 
form.  Usually  striations  manifest  themselves  as  moving 
waves.  There  were  a lot  of  attempts  to  describe  this 
phenomenon  (see,  for  example,  [1]).  However,  in  most 
papers  stratification  was  considered  under  the 
assumption  of  quasineutrality  of  the  positive  column, 
and  the  electric  field  was  introduced  from  some 
additional  consideration  without  taking  into  account  the 
Poisson  equation. 

The  stratification  of  a discharge  is  a collective 
phenomenon  of  mutual  influence  of  electron  and  ion 
plasma  components.  In  the  paper,  electron  kinetics  is 
considered  on  the  basis  of  non-local  Boltzmann 
equation  coupled  with  the  non-stationary  ion  balance 
equation  and  Poisson  equation. 

2.  Model 

A plane  low-pressure  discharge  (p<l  Torr)  is  considered 
in  self-consistent  way.  The  influence  of  edge  effects  is 
neglected  and  the  electric  field  is  assumed  to  be  axial. 
The  non-local  Boltzmann  equation  in  two-term 
approximation  for  the  isotropic  part  of  electron  energy 
distribution  function  (EEDF)  written  in  total  energy- 
space  coordinate  representation  was  used: 
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where  Ng  is  gas  density,  U is  kinetic  energy  of  electron, 
e=U-e0(p(z)  is  total  electron  energy,  <p(z)  is  potential 
distribution  in  the  positive  column  (PC)  of  the  plane 
discharge  between  the  edge  of  the  PC  (z  = 0 cm)  and  the 
anode  (z^O  cm). 

For  ions,  the  non-stationary  continuity  equation  in  the 
drift  approximation  was  considered: 

at  dz 

Vj( z)  is  the  direct  ionization  frequency,  which  in  this 
paper  is  omitted,  |Xj  is  ion  mobility  coefficient. 


The  coupling  of  the  electron  and  ion  components  is 
provided  by  the  Poisson  equation  for  self-consistent 
electric  field 

~fc’  = 4 7R?0  («,-«„),  (3) 

For  the  solution  of  the  system  (1-3),  the  iterative 
numerical  procedure  was  used.  To  obtain  the  first 
approximation  of  distribution  function  /^(z,U),  the 
parabolic  equation  (1)  was  solved  as  an  initial  boundary 
value  problem  from  higher  to  lower  total  energies  in 
some  proposed  electric  field  £°( z).  The  numerical  code 
is  based  on  the  solution  of  tridiagonal  linear  system  of 
discrete  form  of  Eq.  1 , which  is  analogous  to  the  Crank- 
Nikolson  scheme  [2].  Usually  the  constant  field  was 
chosen  as  the  initial  one.  From  the //*(z,U),  the  electron 
density  distribution  n^z)  was  obtained.  With  the  help 
of  Eqs.(2-3),  a new  approximation  of  electric  field  E( z) 
was  found,  and  the  procedure  was  repeated  until 
successive  iterations  of  the  electric  field  were 
converged.  Usually  ten  iterations  were  enough  for  the 
procedure  convergence.  The  final  axial  electric  field 
distribution  E( z)  did  not  depend  on  the  initial 
approximation  £0(z). 


3.  Results 

As  a result  of  numerical  calculations,  spatial  evolution 
of  the  isotropic  part  of  EEDF  was  obtained  (see,  Fig.l). 


Fig.l.  Three-dimensional  representation  of  the  isotropic 
EEDF  starting  from  the  cathode  side  edge  (z=0)  of  PC 
toward  the  anode  (averaged  electric  field  E=4  V/cm, 
p=0.5  Torr). 
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It  was  found  that  the  iterative  process  is  converging 
only  for  the  case  of  moving  striations.  In  plane 
discharge,  the  solution  for  ion  density  and  electric  field 
were  looked  for  in  the  form  E=E(Vjt-z),  with 

the  striation  velocity  ajii£0,  0<a<I. 

The  velocities  of  moving  striation  have  the  value  of 
averaged  ion  drift  velocity,  vj  = ajij£0,  with  parameter 
0<a<l.  For  a<0.25  the  used  iterative  procedure  does 
not  converge.  With  the  increase  of  ot>0.25,  the  peaks  of 
electric  field  in  striations  are  decreased. 

In  Fig.  2,  an  example  of  self-consistent  field  distribution 
is  presented.  It  is  seen  that  electric  field  distribution  has 
non-sinusoidal  peak-like  structure.  The  electron  density 
distribution  is  almost  in  anti-phase  with  respect  to 
electric  field. 
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Fig.2.  Electric  field  and  electron  density  distributions  in 
plane  striated  argon  glow  discharge  (a=0.5). 

It  is  found  that  for  argon  pressures  p>1.5  Torr  striations 
are  damped  due  to  energy  losses  in  elastic  collisions. 
For  lower  pressures,  elastic  losses  result  in  small  decay 
of  striations  to  the  anode  side  of  PC. 
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Fig.3.  Striation  length  dependence  on  reduced  electric 
field  in  plane  glow  discharge. 

In  Fig.3,  the  dependence  of  striation  length  on  reduced 
electric  field  is  shown.  Results  calculated  for  different 
argon  pressures  and  electric  fields  fall  into  the  unified 


curve.  For  given  E/p,  the  striation  length  is  inversely 
proportional  to  gas  pressure,  L~I/p.  The  length  L 
slightly  exceeds  the  minimal  length  necessary  for  an 
electron  to  attain  the  first  excitation  threshold 
L>Ui/e0E0. 

In  Fig.  4,  the  distribution  of  volume  charge  is  shown  for 
averaged  field  £=4  V/cm  and  p = 0.5  Torr.  It  is  seen 
that  double  layers  are  formed  in  stratified  PC.  Relative 
deviation  from  the  quasineutrality  has  the  order  of  10*2- 
10*3. 
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Fig.4.  Volume  charge  distribution  in  stratified  PC. 

4.  Conclusions 

It  this  paper,  a self-consistent  model  of  moving  striation 
in  argon  is  presented.  Non-local  Boltzmann  equation, 
non-stationary  ion  balance  equation,  and  Poisson 
equation  were  considered  simultaneously  with  the  help 
of  iterative  procedure. 

It  was  found  that  for  low-pressure  gas  conditions  the 
stratified  regimes  cannot  be  realized  for  the  case  of 
standing  striations.  However,  consideration  of  moving 
striations  permits  us  to  receive  a converged  solution  of 
the  problem.  Only  the  moving  striations  with  the 
velocities  directed  from  the  anode  to  the  cathode  adjust 
the  phase  shifts  between  electric  field  and  electron 
density  distribution.  The  resonant  amplification  of  peaks 
of  distributions  takes  place. 

The  presented  kinetic  model  is  applied  for  low-pressure 
inert  gases.  At  intermediate  pressures,  the  stepwise 
ionization  from  metastable  states  should  be  taking  into 
account.  In  this  case,  system  of  Eqs.(l-3)  has  to  be 
supplemented  with  corresponding  balance  equation  for 
matastable  particles 
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